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We study the statistical properties of Ising spin chains with finite (although ar- 
bitrary large) range of interaction between the elements. We examine mesoscopic 
subsystems (fragments of an Ising chain) with the lengths comparable with the in- 
teraction range. The equivalence of the Ising chains and the multi-step Markov 
sequences is used for calculating different non-additive statistical quantities of a 
chain and its fragments. In particular, we study the variance of fluctuating magne- 
tization of fragments, magnetization of the chain in the external magnetic field, etc. 
Asymptotical expressions for the non-additive energy and entropy of the mesoscopic 
fragments are derived in the limiting cases of weak and strong interactions. 
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I. INTRODUCTION 

One of the basic postulates of the traditional statistical physics is the assumption that 
the range of particle interaction is small compared to the system size. If this condition does 
not hold the internal and free energy, entropy, and other physical quantities are no longer 
additive. Due to this fact, there are no straightforward definitions of the temperature, 
entropy, etc. The statistics of the systems with long-range interaction are not Gibbsian; 
Boltzman's relationship between the entropy and the statistical weight is no longer valid; 
the fluctuations of sums of random quantities do not obey the Gaussian statistics. 

Systems with long-range interaction are common in physics Amongst them are: 
gravitating systems electrically charged systems {3)], magnets with dipolar interactions {4], 
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etc. An important physical example of such class of systems is the Ising spin chain with 
long-range interaction (see, e.g., [5|). 

The long-range interaction leads to the correlations between the elements of the system, 
which is the reason for the appearance of non-Gaussian statistics. Such non-Gaussian distri- 
butions (induced by the long-range correlations) are found in various fields of science, e.g., 
in linguistics (distribution of "words" in a given text , biology (the distribution of 

nucleotides in DNA sequences [9|), computer science (forward error correction codes 10]). 
economics (the distributions in financial markets 111), sociology, physiology, seismology. 

One of the most interesting phenomena in the long-range correlated systems is the phe- 
nomenon of phase transition. The existence of phase transitions is of great importance for 

Tin n 

different sciences, for example, biology in the problem of DNA melting |12l. Il3|. Dyson [14j . 



Ruelle 



15], Thouless [16|] found different conditions for existence of the phase transition in 



the Ising chain with the infinite scale of interaction. 

The usual objects for studying the effects of long-range correlations are the non-extensive 
thermodynamic systems where both interaction range and the system length are macro- 
scopically large. However, the mesoscopic long-range correlated systems, which are not large 
enough to consider the thermodynamic limit, have not been studied thoroughly. The statis- 
tical properties of such systems are very similar to the properties of the infinite non- extensive 
thermodynamic systems. In the present paper, we consider the finite-size subsystems (frag- 
ments) of an infinite Ising chain as examples of mesoscopic correlated systems. We study 
the Ising chain with the Hamiltonian, 

H = - e(] ~i)siSj - H^Si, (1) 

0<j—i<N i 

where the spin variable s« takes on two values, —1 and 1, e(r) > is the exchange integral of 
the ferromagnetic coupling, and H is the external magnetic field. Below, for the simplicity, 
we assume e(r) = for r > N. The range N of spin interaction is finite although arbitrary. 
The chain is assumed to be in the thermodynamical equilibrium with the Gibbs thermostat 
of temperature T (unless otherwise mentioned). Such a system is obviously extensive one, 
and the usual thermodynamics can be applied to the system as a whole. There is no 
spontaneous magnetization in such a chain due to the finiteness of the particle interaction 
range. Nevertheless, its mesoscopic subsystems are of interest because of their non-additive 
statistical properties. The aim of our work is to calculate different statistical quantities for 
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an mesoscopic subsystem of length L, that can be of the order of the interaction length N. 

Thus, we study an Ising chain with the finite range of interaction, but the lengths of 
the subsystems are also finite. We show that energy and entropy of the fragments of the 
Ising chain are non-additive quantities. We calculate the root-mean-square (RMS) value of 
fluctuation of the magnetization which does not scale as \[L (the number of spins in a frag- 
ment), unlike the case of short-range interactions. We suggest a way to calculate statistical 
quantities for the finite subsystems in contact with both the external Gibbs thermal bath of 
temperature T and the rest of the Ising chain playing a role of additional thermostat. 

For mesoscopic systems, an important problem is to introduce the appropriate statistical 
ensemble correctly, because different statistical ensembles are not equivalent for such systems 



(see, e.g., Ref. Il8l). These canonical ensembles are equivalent in the thermodynamical limit 
only. In the present paper, we deal with an ensemble of the finite subsystems that are the 
parts of an infinite Ising chain. 

We calculate the statistical properties of the Ising chains using their equivalence with the 
binary N-step Markov chains. The multi-step Markov chain is the sequence where symbols, 
say ±1, are generated by means of the following procedure. One starts from arbitrary N 
symbols, than the subsequent symbols are generated with the probability that depends on 
the values of the preceding N symbols only and is independent of the farther symbols, 

P(si = s\T.) = P( Si = s\T- N ). (2) 



Here T i L is a set of L sequential symbols (si-L, s 



of such chains were studied in detail in Refs. 
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Vi). The statistical properties 
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Statistical equivalence of the Ising and Markov chains is widely known (see, e.g., Refs. 
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25l ). However, these References define Markov chains not via the "one-sided" conditional 
probability function Eq. (j2J), depending on the values of the preceding symbols, but rather 
via the "two-sided" one. The definition of the chain with one-sided conditional probability 
function is more convenient for numerical simulations of the chain as well as for obtaining 
analytic results. Reference 



26l demonstrated the equivalence of the two mentioned views on 
Markov chains. 

The paper is organized as follows. In the second Section we show how to find the Markov 
chain corresponding to the Ising chain under consideration. Using this correspondence rule 
we calculate the RMS value of the fluctuating part of magnetization for the chain fragments 
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of the length L in the absence of magnetic field. We also find magnetization of the chain in 
the presence of the magnetic field. The third Section introduces the method for calculating 
the statistical quantities of the fragments of Ising chain, in particular, the internal and total 
energy, and entropy. All these quantities are non-additive, i.e., they do not scale as the 
fragment length L. We obtain analytical results for the two limiting cases of low and high 
temperatures. 



II. MAGNETIC PROPERTIES OF ISING CHAIN 

First, we find a one-sided conditional probability function of an iV-step Markov chain, 
Eq. (J2D, that is statistically equivalent to the Ising chain with Hamiltonian ([T|). Then 
we use this property to find the magnetization of the Ising chain for the case of weak 
interactions (compared to the temperature). Finally, we study another limiting case of 
strong interactions. 



A. Equivalence of Markov and Ising chains 



The conditional probability of a spin to have a definite value, say 1, given that the values 
of all other spins in the chain are fixed, is independent of the values of the spins separated 
by the distances larger than N (see, e.g., Refs. 
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25l ) and reads 



P( Si =l\T-T^ ) 



P( Si =l\Tr N ,T, 



Here T^ L is the L-tuple, (sf+i, s i+2 , ■ ■ ■ , s i+ L). The probability function in Eq. 
by the Gibbs formula, and for the Hamiltonian ([TJ) it has the following form: 



JV 



P( Si = l\T- m T+) 



, 2H ^2e(r,. 

1 + exp ( - — - 2^ -J^{Si-r + Si+r)) 



r=l 



(3) 



is given 



(4) 



Expression (j4j) for the two-sided conditional probability function can be derived by the 



method proposed in Section III. It can be shown that the Metropolis scheme 
the same expression (T4j) for the conditional probability. 



271 ] also yields 



Reference ]26| analytically proved that the chain defined by a two-sided conditional proba- 
bility function (]3]) is equivalent to the iV-step Markov chain defined by Eq. (T5]). The relation 
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between the two-sided and one-sided Markov conditional probability functions reads 

N 

EI P( S i+r | T i+r ^ N ) 

r=0 

P( Sl = s I Tr m T+ N ) = — ^ . (5) 

II P(s i+r | Tr ) + n P(s i+r I T- r ^ N ) 

r=0 r=0 

Sj = l S< = — 1 

Unfortunately, this relation is quite cumbersome and can hardly be applied to the an- 
alytical study of Ising chains in the general case. Nevertheless, it allows one to calculate 
numerically the one-sided conditional probability function for the iV-step Markov chain. 
Thus, to attain equilibrium state of the spin chain, we can generate the Markov chain 
according to the algorithm described in the Introduction, instead of using the Metropolis 
scheme. 



B. Magnetization of the Ising chain with weak interaction 

In order to use Eqs. (jlj) and ([5]) for the analytical study of Ising chains we consider the 
case of weak interaction compared to the temperature, 

A' 



£e(r)<T. (6) 



r=l 

N 

Expanding Eq. (jlj) over the small parameter Y2 £ ( r )/T and using Eq. ([5]) we obtain the 

r=l 

one-sided conditional probability function. In the first approximation, it takes the additive 
form, 

N 

P(s t = 1 | T-J = ^ + J2 F (r)(si-r - s). (7) 

r=l 

Here the function F(r) and s (the average value of Sj) are determined as follows, 

F(r) = ^cosh- 2 ^, S -=(l + 4f>(r))tanh^. (8) 

r=l 

The Markov chain with additive conditional probability function, Eq. ([7]), is referred to as 
the additive Markov chain. The function F(r) is referred to as the memory function 
Equation ([8]) reflects the obvious fact that an increase of the temperature leads to an increase 
of disorder and to decrease of correlations in the chain. It should be noted that the interaction 
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range N in Eqs. (J7|) and (jSJ) can be taken as oo, if the series ^ r converges and the 

value of P in Eq. ([7j) satisfies the inequality < P(sj = 1 | T^) < 1. 

Thus, the Ising chain in the equilibrium is equivalent to the additive Markov chain. For 
the case of weak interaction between the spins, the corresponding memory function F(r) 
is proportional to the energy e(r). The memory function of the additive Markov chain is 
related to its pair correlation function K(r), 



K(r) = SiS i+r - s , (9) 



by the recurrence equation 28 



K(r) = ^2F(r')K(r-r'), r = 1, 2, . . . . (10) 

r'=l 

The overline ... in Eq. (J9j) denotes statistical averaging. In the limiting case of small F(r), 
Eq. (TTTJT) yields approximately, 

K(r)&2F(r), r = l,...,N. (11) 



1. Fluctuations of magnetization in the absence of magnetic field 

We now consider an Ising chain with weak but long-range interaction in the absence of 
the magnetic field. As mentioned above, the averaged magnetization equals to zero in this 
case. However, the root-mean-square (RMS) value of the magnetization of a chain segment 
with length L, 

/ \ 1/2 

M RMS (£) = [{si+i + s i+2 + ■■■ + s i+L ) 2 j , (12) 
is not zero and is determined by the correlation function of the chain, 

1/2 

M RMS (L) = [L + 2 - r)AT(r)) . (13) 

r=l 

Using Eqs. (JHJ) and (jlip . we derive the relation between the RMS of magnetization and 
the interaction energy, 

1 L_1 

M RMS (L) = VI + _ J^iL - i)e{i). (14) 

* i=l 

Here we used that e(r) = for r > N, thus the range TV of interaction does not appear 
explicitly. 
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The interaction between spins leads to the deviation of the RMS of magnetization from 
the value a/L occurring in the non-interacting chain. If this deviation is small, as we assume 
here, the distribution function of magnetization is close to the Gaussian form with the 
variance M^ MS (L). 



2. Magnetization in the external magnetic field 

If the magnetic field H is applied, the averaged value \x of magnetization (per spin) is 
no longer zero. In the assumption of weak interactions, Eq. (J6]), fi is equal to the averaged 
value of the spin s, Eq. (JSJ) , 

N 

V = s= (^l + -cosh- 2 — ^e(r)Jtanh— . (15) 

r=l 

In the nearest-neighbor approximation (N = 1), this equation coincides with the high- 
temperature asymptotics of the well-known expression for the magnetization, 



exp(-4e/T)\ 



-1/2 



^^ H { 1+ ^?mM) ■ (16) 

Here e = 5(1) is the single exchange integral for the nearest neighbors. 

The results of numerical simulations for the Ising chain with weak interaction are given 
in the left panel of Fig. [TJ Here we present the dependence of the magnetization on the 
value of magnetic field. The inset shows the dependence of the RMS of magnetization on 
the segment length L in the absence of the magnetic field. The numerical points are very 
close to the corresponding analytical curves (solid lines) determined by Eqs. ( fl5l) and ( fl4l) . 



C. Magnetic properties of the Ising chain in the low-temperature limit 

In the previous subsection, we have used the statistical equivalence of the Markov and 
Ising chains for the calculation of the magnetization in the case of high temperatures. Here 
we examine the opposite limiting case of low temperatures. The corresponding inequality 
imposed on the temperature T, the interaction-energy function e(r), and magnetic field H 
will be given below (see Eq. ( fl9i) ). In this case, Eq. (J4]) defines a non- additive Markov chain 



fsee Ref. 
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for details). Since finding the correlation function for such chains brings about 
considerable technical difficulties, we examine the required statistical properties of the Ising 
chain directly from Eq. fflh. 
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FIG. 1: (Color online) The magnetization fj, of the chain per spin versus the magnetic field H. 
Left panel (T=10): solid and dashed curves correspond to the high-temperature, Eq. (|15fl . and 
low-temperature, Eq. (|22p . limits, respectively. Right panel (T=0.6): solid curve corresponds 
to low-temperature, Eq. (|22p . limit. The RMS magnetization in the absence of magnetic field 
versus length L of a segment is presented in the insets. The dots correspond to the results of 
the numerical simulations. The function e(r) is step-like, e(r) = 1/100 for r = 1,2, ... , 100 and 
e(r) = for r > 100. 

1. Statistical properties of strongly correlated Ising spin chain 

In the case of low temperatures, the conditional probability determined by Eq. (jlj) is 
close either to zero or one. Indeed, due to strong correlations, the chain consists of large 
domains predominantly occupied by the same symbols. It is evident that, once strong enough 
magnetic field is applied, the chain will mostly consist of the same symbols. 

For definiteness, we consider the case when almost all symbols in a chain are 1 (the case 
of if > 0). In other words, the probability to find the spin —1 at a given place of the chain 
is small. For further convenience, it is suitable to introduce the parameter 5, 

N 

5 = exp(-4ei/T) < 1, e 1 = ^e(r) + if/2. (17) 

r=l 

According to Eq. (jlj), this parameter is a ratio of the probability to find the (2N + l)-tuple 
with only one spin —1 in the center to the probability of occurring the (2N + l)-tuple with 
all spins 1. Due to the correlations, the probability to have the (2N + l)-tuple with one 
more spin —1 at distance k from the central spin —1 differs from 5 2 and is equal to 5 2 ak 
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with 

a k = exp(4e(£;)/T) > 1. (18) 

Further, due to Eq. (J4]), three spins —1 (at successive distances k and I) occur even less 
frequently, with the probability 5 3 <Jk<Ji<Ji+k- The corresponding probabilities for four and 
more spins —1 can be found in a similar way. 

We use the following assumptions for the parameters of the system: 

N 

Sa^af < (1 +exp{-2H/T)y 2N . (19) 
i=i 

This inequality states that the probabilities to have a (2N + l)-tuple with three or more 
spins —1 is negligible compared to that with just two spins —1. Inequality (Tl9|) holds only 
for the case of monotonously decreasing functions e(r). In general case, we should replace 
G\ of in Eq. (fT9l by the maximum value of a k aiai+k- 

Inequality ( fl9i) can be considerably simplified under some suppositions. For example, if 
function e(r) is rapidly decreasing (e(l) ^> e(2) ^> . . . ) and H > 4e(l) we obtain 

T<#-4e(l). (20) 

For chains with approximately step-like function e(r) (e(r) » 5(1) for r < N and e(r) = 
for r > N), we have 

T Ne(l) + H/ h\N. (21) 

Thus, Eq. ( Fl9l) means that the temperature T should be sufficiently small compared with 
H or e(i). 

It should be noted that inequalities ffl9l) and ([6]) can be satisfied simultaneously. This is 
the case when ^ sir) <T< H/N. Then, the results of subsection III Bl coincide with those 
of this subsection (Eq. ( fl5l) transforms to Eq. (J22j)). 

2. Magnetization of the Ising chain 

Using Eqs. (j4j), (fl71) . and (TTBl) . we derive average magnetization fi per spin in the presence 
of the magnetic field H > 0, 

N 

f i = s = l- 25 + 25 2 - 45 2 J2( ( rk-l)- (22) 

k=i 
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For strong enough magnetic fields, exp(— H/T) <C 1, this expression is also valid in the high- 
temperature limit Eq. ([6]). Moreover, for chains with nearest-neighbor interaction (N = 1), 
it coincides with the general formula ( fT6l) in the main and first approximations with respect 
to the parameter exp(— (H — 4e)/T) <C 1. 

If the magnetic field is not applied, the averaged magnetization is equal to zero. However, 
the RMS of magnetization, determined by the pair correlation function K(r), is non-zero. 
In the case of low temperatures, K(r) is close to one. For r ^ N, 

K(r) « 1 - 45. (23) 

Therefore, the RMS of magnetization of segments with length L <C L is equal to L in the 
main approximation with respect to 8, contrary to the square-root dependence \J~L valid 
in the high-temperature limit. Here Lq ^> N/8 is the characteristic length of the domains 
consisting of the same symbols, 1 or —1. Using Eqs. (TT31 and (1231 . we obtain, 



M BM3 (L)=L-2{L-1)6, (24) 

The results of numerical study of strongly correlated chains are shown in Fig. [U As in 
the case of weak interaction, we present calculations of the RMS of magnetization for the 
segment of length L in the absence of the magnetic field (inset in the right panel) and the 
averaged magnetization as the function of the magnetic field H (solid curve in the left panel 
and dashed curve in the right panel). All the results of numerical simulations are close to the 
theoretical predictions. As mentioned above, two asymptotics Eqs. (Tl5l) and (122j) coincide 
for the strong enough magnetic fields at Yl e(r) <C T. This is clearly seen in the left panel. 

Thus, we have studied the magnetization and its RMS for the Ising chains in the two lim- 
iting cases of high and low temperatures. These quantities are determined by the statistical 
properties of the chains, particularly, by the correlation function. Not only the magneti- 
zation, but all the other statistical properties of the Ising chain can be studied using the 



equivalence of the Ising and Markov chains (see Ref. l2ll ). In the next section, we propose a 
general way of calculating the statistical quantities of the Ising chains with long-range spin 
interaction. 
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III. TOTAL AND INTERNAL ENERGY AND ENTROPY OF ISING CHAIN 

SEGMENTS 

Here we present an approach to calculate any statistical quantity (energy, entropy, etc.) 
for a subsystem of arbitrary length L (a segment of L sequential spins in the Ising chain). 
This mesoscopical subsystem, denoted by M, interacts with the thermostat and with the 
rest of the chain. We denote also two "border" subsystems of the length N by B. All the 
other spins, except for M and B, are denoted by R: 



■ ■ ■ Si—N S i-N+1 ■ ■ ■ Si Sj+i . . . Si+L s i+L+l ■ ■ ■ S i+L+N S i+L+N+1 

-V ' N v 



R B M B R 

The total energy of the chain is not equal to the sum of the internal energies of subsystems 
M and (B + R). So, the long-range interaction between M and (B + R) makes it impossible 
to use directly the standard methods of statistical physics for calculation the statistical 
quantities for segment M. 

In order to find the condition of equilibrium, we introduce the statistical ensemble of 
the chains. We choose the subensemble of chains with fixed border regions B (but with 
varying M and R). For each of 2 2N such subensembles, the subsystem B plays the role of 
a separating wall. If the ensemble is in the equilibrium state, the total energy of (M + R) 
(that includes the sum of internal energies of M and R, and energy of their interaction with 
B) is constant for every subensemble. Thus, we can use the equilibrium condition between 
M and R: 

d\nW M (E M \B) = dlnW R (E R \B) = 1 

dE M dE R T ' 1 j 

where Wm(E m \B) and Wr(E r \B) are the statistical weights of subsystems M and R with to- 
tal energies Em arid Er, with the borders B fixed. We refer to statistical weights Wm{Em\B) 
and Wr{E r \B) as the conditional statistical weights. It should be emphasized that these 
weights are not mesoscopical quantities because of their dependence on the microscopic 
states of borders B. If the system is in the thermal contact with the Gibbsian thermostat, 
the condition of equilibrium between the thermostat and segment M sets the temperature 
of M equal to T. This temperature is obviously the same for every subensemble with fixed 
B. Thus, the averaged temperature of the subsystem M is also T. 
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Note that the external thermostat is not necessary for establishing the temperature T of 
the mesoscopical segment M. The role of a thermostat can be played by the subsystem R. 
Indeed, the temperature determined by Eq. (1251) does not depend on values of spins in B 
(because the statistics of infinite subsystem R can not depend on the microstate of finite 
subsystem B). This means that, for all the subensembles, the temperature of M is the same 
even in the absence of the external thermostat. In this case, the temperature in the chain 
is determined by its initial state (that can be nonequilibrium one). 

It is very important that the distribution function of the segment M over different mi- 
crostates within a subensemble with fixed borders B is Gibbsian (though the segment M 
itself is not a Gibbsian system). Therefore, within every subensemble, we can introduce a 
conditional statistical quantity Qm(\B) for the segment M. The actual quantity Qm is the 
conditional one, averaged over the subensembles with different borders B: 



With this method for calculating the statistical quantities, we do not need to find the 
distribution function for the segment M over different microstates in the whole ensemble. 
Note also that Eq. (1261) can be considered in the thermodynamical limit L — > oo. However, 
we focus our attention on the mesoscopical segments M with finite length L. 

The conditional entropy can be introduced as the logarithm of the conditional statistical 
weight: S M (E M \B) = \nW M (E M \B). Equality dS M (E M \B) = dE M {T\B)/T is fulfilled for 
the conditional quantities Sm{\B) and Em{\B). Meanwhile, such a relation is not valid for 
the averaged entropy S M {E M ) = (S M (E M \B)) B and total energy E M (T) = (E M (T\B)) B . 

Note that the presented method for calculating the statistical quantities is rather general 
and can be applied, e.g., to the internal energy of a segment M that can be measured 
experimentally, or to the probability for some spin to take on the definite value under 
the condition of fixed environment. Considering one spin Sj as a subsystem M and the 
conditional probability P as a quantity Q in the above-mentioned method, one arrives at 



Now we apply our approach to calculation of the total and internal energy and entropy 
of the segments of length L for the limiting cases of high and low temperatures. 




(26) 



Eq. ®. 
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A. High-temperature limit 

The total energy of a segment M with length L for prescribed configuration of M and B 
is expressed via products of pairs of spins, 

i+L i+L j+N i+L 

E(L\M,B) = -H s j~ E E z(\k-j\)s kSj + Yl <k-j)s k s r (27) 

j=i+l j=i+l k=j—N j,k=i+l 

k^j j<k 

The actual total energy is the conditional energy averaged over the different configurations of 
M and B, E(L, T) = E(L\M, B). Along with the total energy E(L, T), the segment M can 
be characterized by its internal energy E m (L,T) = E m (L\M, B) specifying the interaction 
between the spins in M only, 

i+L i+L 

E' m (L\M,B) = -H Sj~ E e(k-j)s kSj . (28) 

j=i+l j,k=i+l 
j<k 

This quantity can be measured experimentally and thus is of interest for study. 

The second term in Eq. (1271) accounts the interaction of each spin of the segment M 
with 2N surrounding spins. The internal interaction energy is accounted twice in this term. 
Thus, the third term in Eq. (1271) subtracts the corresponding extra summands. 

In the high-temperature limit, the total and internal energy can be calculated via the 
pair correlation function K(r) with r ^ A, without the need for the conditional energies. 
Using Eqs. (jHJ) and (1TT1) . we arrive at 

L-l 

E(L,T) = -{e 2 L - ~ r)e\r))/T, (29) 

r=l 

L-l 

E in (L, T) = - ( £(L - r)e\r) + H 2 L) /T, (30) 

with 

N 

e 2 = 2^e 2 (r) + i/ 2 , #<T. (31) 

r=l 

If the segment length L is much greater than the memory length A this expression yields 
additive energies, E(L,T) oc L. For the opposite limiting case, L A, we get: 

E(L, T) w -e 2 L/T + e 2 (l)L 2 /2T, E in (L, T) w -5 2 (1)L 2 /2T - // 2 L/T. (32) 

Here we suppose that e(i) ~ e(l) for z < L. 
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FIG. 2: (Color online) The dependences of specific non-additive energy E(L)/L and E m (L)/L on 
the length L of a segment for the step-like function e(r) = 1/100, r < 100. The specific non- 
additive entropy S(L)/L versus length L is presented in the inset. Other parameters are: T = 10, 



The non-additive energy is expressed in terms of the pair correlation function only. This 
is not correct for other statistical quantities, e.g., the entropy of the segment M. Formally, 
in order to find the entropy, one should calculate all the conditional entropies by integrating 
the equation dS(L,T\B) = dE(L,T\B)/T, and averaging the result over all realizations 
of the borders B. However, at high temperatures, to a first order in the small parameter 
e\/T, we can change the order of these operations and calculate the averaged entropy by 
integrating the formula dS(L,T) = dE(L,T)/T written for the averaged energy. A con- 
stant of integration is determined from the condition of complete randomization at high 
temperatures, S(L,T — > oo) — > In (2 L ). Thus, we obtain 



with E(L, T) given by Eq. fl29|) . 

Expressions (|29|) - (!33|) describe the non-additive dependence of energy and entropy of a 
segment M on its length L. Note that the energy is non-additive in the main approximation 
in the parameter ei/T, while the non-additivity of the entropy appears in the first-order 
correction only. The dependences of the non-additive energy and entropy on L are presented 
in Fig. [2] for the step-like function e(r). 

Knowing the energy and entropy we can find some other statistical quantities. For exam- 
ple, at high temperatures the heat capacity can be determined in a similar way as for the 



H = 0.1. 



S(L, T) = L In 2 + E(L, T)/2T 



(33) 
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entropy. One can use the classical formula C(L,T) = TdS(L,T) jdT with averaged entropy 
given in Eq. fl33|) and obtain the equation 

L-l 

C(L,T) = -E{L,T)/T= [e 2 L-J2(L-r)e 2 (r))/T 2 . (34) 

r=l 

This formula is valid for the case of high temperatures in the main approximation only. In 
the general case, C(L,T) is not determined by equation C(L,T) = TdS(L,T)/dT. This 
relation holds for the conditional quantities only. 

For the case of high temperatures, we calculated the averaged non-additive energy and en- 
tropy without using the conditional ones. In the opposite limiting case of low temperatures, 
it is necessary to calculate the conditional quantities as well. 



B. Low-temperature limit 

In this subsection, we assume that inequality Eq. (TT9]) is fulfilled. Let us find the averaged 
values of the total and internal energy and entropy of a segment M consisting of L spins. 
These are close to their minimal values taken at T = when all spins in segment M (as 
well as in borders B) are oriented along the magnetic field H . For positive magnetic fields, 
H > 0, we have 

N N 

E min (L) = -(ex + H/2)L - J>e(r), E^(L) = -{e x + H/2)L + £Ve(r). (35) 

r=l r=l 

For simplicity, we consider here and below the segments with L ^ N. 

According to the Nernst theorem, the entropy in the zeroth approximation is S(L, T = 
0) = 0. Now we seek the first non-additive corrections to E, E m , and S. Following the 
proposed approach, we have to find, as a first step, the conditional total energy and entropy 
for all the configurations of borders B: 

i) The most probable configuration of borders B is Bq when they consist of spins +1 only. 
The corresponding probability P(B ) is approximately equal to 1 — 2N5. The conditional 



16 



total energy and entropy for this case can be found by virtue of Eqs. (fT7|) . ({TBI , and (I27|) : 

L 

E{L, T\B ) = E min {L) + 4<feiL(l - L5) + 5 2 ^ (8ei - 4e(A;)) (L - k)a k , (36) 



fc=i 

T 

S(L,T\B ) = J ^dE(L,T\B ) 
o 

= « jfi+l - ^gfl±l + ^ ± * ~ *f + \ L - tK ,37, 

fc=l 

ii) Much less probable configuration of the borders is B\ k with only one symbol —1 at 
distance k from the edge of segment M. The corresponding probability does not depend on 
k and is approximately equal to 5, P(Bi^) ~ 5. For this configuration, the conditional total 
energy and entropy are: 

L L 

E{L, T\B ltk ) = E min (L) +2j2<k + j-l) + 8j2 ( 4e i ~ + j ~ (38) 

j=i j=i 

S(L, TIB,,) = j \iE(L, T\B ht ) = * £ * ~ + ' ~ » + T ^ 







3=1 



iii) Probabilities of all other configurations of the borders are negligible with respect to 
P(B\^)- There is no need to find the corresponding conditional energies and entropies if we 
are interested in the first non-additive corrections to these quantities. 

Averaging the conditional quantities over border configurations B and B\ k and keeping 
the non-additive temperature dependent terms in the main approximations in 5, we derive 
the actual quantities E(L,T) and S(L,T): 

N 

E(L, T) = E min (L) + re{rj) , (39) 



r=l 

N 



S(L, T) = LS [£(1 - 26) + 1 + + S £ *--f> + V - 1)] (40) 

fc=l 

AT . AT 

^[E^-^-rE^H- (4i) 



fe=l k=l 

Similar calculations for the internal energy give the following expression: 

N 



E%{L,T) = ^(L) + 4tf(eiL - J>e(r)). (42) 

r=l 
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FIG. 3: (Color online) The specific non-additive energies E, E m and entropy S versus the length 
L of a segment for step-like function e(r) = 1/20, r < 20. Symbol A denotes the non-additive 
terms in Eqs. (i39j) . (f4T1) and (f42j) . Other parameters are: T = 3, = 6. 

One can see that the non-additive total energy and entropy, determined by Eqs. ( 1391) 
and fid]) , do not satisfy the relation dE(L,T) = TdS(L,T), unlike the conditional total 
energy and entropy. The dependences of the non-additive energy and entropy on L are 
presented in Fig. [3] for a step-like function e(r). 

As mentioned above, in the absence of an external thermostat, the temperature T of 
a chain is controlled by its initial state. In this case, the temperature can be expressed 
via the total energy of the chain. In particular, the positive initial energy of the chain 
provides negative temperature (see Fig. [4]). If the initial energy is close to its maximum, the 
temperature is negative and close to zero and the system is almost "frozen" . The lesser the 
positive energy the lower the negative temperature. If the initial energy tends to zero, the 
chain becomes more chaotic and its temperature tends to the minus infinity. 

The situation is principally different for the chains in contact with an external thermostat. 
The states of chain with negative temperatures become unstable and, therefore, do not exist. 
If the initial energy of the chain is negative and close to zero, the temperature is positive 
and high. At the same time, it tends to zero for the energy close to its minimum. Both with 
and without the external thermostat, all the statistical properties of the chain with positive 
temperature are the same. 
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FIG. 4: (Color online) The dependence of the entropy S on the total energy E for the step-like 
function e(r) = 1/20, r < 20. Magnetic field H = 100 is so strong that the high-temperature (red 
dashed curve), Eqs. ([29]) . ([33]) . and low-temperature (blue solid curve), Eqs. ([39]) . (|4Tj) . asymtotes 
coincide. The corresponding values of the temperature are shown on the right scale. 



IV. CONCLUSION 



Thus, we have studied the statistical properties of mesoscopic segments of Ising spin 
chains with finite but arbitrary long interaction ranges. The equivalence of the Ising and 
iV-step Markov chains was used for calculating the averaged statistical quantities. In partic- 
ular, the averaged magnetization in the presence of the external magnetic field and its RMS 
were calculated in the two limiting cases of strong and weak interactions. Correlations in 
the chain result in the non-additive behavior of these quantities as the segment length L in- 
creases. We show that the statistical quantities of the chain can be obtained by averaging the 
corresponding conditional quantities. The explicit expressions for the non-additive energy, 
internal energy, and entropy are derived in the limiting cases of high and low temperatures 
comparing to the energy of spin interaction. At high temperatures, the equilibrium Ising 
chain of spin turns out to be equivalent to the additive multi-step Markov chain. 
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